NASA Technical Memorandum 104143 


4/ocJO 



UNSTRUCTURED-GRID METHODS DEVELOPMENT FOR 
UNSTEADY AERODYNAMIC AND AEROELASTIC ANALYSES 


a i^7«S!T; ,e c !“ 01A 


N9l-3208<* 


UncTas 

l ’3/02 0046090 


JOHN T. BATINA 
ELIZABETH M. LEE 
WILLIAM L. KLEB 
RUSS D. RAUSCH 


SEPTEMBER 1991 


NASA 

National Aeronautics and 
Space Administration 

Langley Research Center 

1 tampion. Virginia 23665 


f 



IJNSTRIJCTURED-GRID METHODS DEVELOPMENT FOR 
UNSTEADY AERODYNAMIC AND AEROELASTIC ANALYSES 

John T, Baltna 
Elizabeth M. Lee 
William L. Kleb 
NASA Langley Research Center 
Hampton, Virginia 23665-5225 USA 

and 

Russ D. Rausch 
Purdue University 

West Lafayette, Indiana 47907 USA 


Summary 

The current status of unstructured-grid methods develop 
ment in the Unsteady Aerodynamics Branch at NASA Langley 
Research Center is described. These methods are being devel- 
oped for unsteady aerodynamic and aeroelastic analyses. The 
paper first highlights (he flow solvers that have been developed 
for the solution of the unsteady Euler equations and then gives 
selected results which demonstrate various features of the ca- 
pability. The results demonstrate two and three-dimensional 
applications for both steady and unsteady flows. Comparisons 
are also made with solutions obtained using a structured grid 
code and with experimental data to determine the accuracy of 
the unstructured grid methodology. These comparisons show 
good agreement which thus verifies the accuracy. 


Introduction 

Considerable progress has been made over the past two 
decades on developing computational fluid dynamics (CFD) 
methods for aerodynamic analysis. 1,2 Recent work in CFD has 
focused primarily on developing algorithms for the solution of 
the Fuler and Navier-Stokes equations. For unsteady aerody- 
namic and aeroelastic analysis, these methods generally require 
that the mesh move to conform to the instantaneous position 
of the moving or deforming body under consideration. Many 
of the methods that arc currently being developed assume that 
the mesh moves rigidly or that the mesh shears as the body 
deforms. These assumptions consequently limit the applicabil- 
ity of the procedures to rigid-body motions or small-amplitude 
deformations. Furthermore, these methods of solution typically 
assume that the computational grid has an underlying geometri- 
cal structure. As an alternative, algorithms have been developed 
recently which make use of unstructured grids. 3-19 In two di- 
mensions these grids are typically made up of triangles and in 
three dimensions they consist of an assemblage of tetrahedra 

The unstructured grid methods have several distinct advan- 
tages over structured grid methods which make them attractive 
for unsteady aerodynamic and aeroelastic analyses. For exam- 
ple, the primary advantage of the unstructured grid methodology 
is die ability to easily model very complicated three-dimensional 
geometries such as the F/A-1K aircraft shown in Fig. I. 9 The 


aircraft was modeled by including the wings with leading edge 
extension, horizontal and vertical tails, as well as the canopy 
and the fuselage. The modeling also includes engine inlets and 
nozzles to simulate engine power effects. With a structured 
grid, it is extremely difficult to achieve this level of geometrical 
complexity. A second advantage is that the methodology allows 
for a general way to move the mesh to treat realistic motions 
and structural deformations of complete aircraft configurations. 
An example of the deforming surface grid for a transport-type 
configuration undergoing a complete-vehicle bending motion is 
depicted in Fig, 2. The deforming grid capability does not 
involve any assumptions which limit applications to small de- 
formations, such as the simple grid shearings done in some 
structured grid codes. A third advantage is that it enables in a 
natural way for adaptive mesh refinement to predict more accu- 
rately the physics for the flow. For example, shown in Fig. 3 
is a conical vortex-dominated flow solution for a flat plate delta 
wing at a supersonic freestream Mach number. 13 The solution 
was obtained by adapting the original coarse mesh three times 
to the instantaneous flow. The final result is a highly accurate 
solution of the Euler equations, produced by using an order of 
magnitude fewer grid points than if a globally fine mesh was 
used. Similar to spatial adaption, temporal adaption may also 
be employed with unstructured grids for unsteady problems to 
resolve more accurately and efficiently the physics of the flow in 
time, 18 Temporal adaption can be thought of as time accurate lo- 
cal time-stepping where smaller time steps are used in grid cells 
where the temporal gradients are large and larger time steps are 
used in cells where the gradients are small. Time accuracy is 
maintained by bringing all grid cells to the same time level as 
determined by the largest step size in the problem. 

The purpose of the paper is to describe the current status 
of unstructured-grid methods development within the Unsteady 
Aerodynamics Branch at NASA Langley Research Center. 10 " 19 
The paper first highlights the flow solvers that have been de- 
veloped for solution of the time-dependent Euler equations and 
then gives selected results which demonstrate various features 
of the capability. The flow solvers that are described are either 
of the central-differencc-type with explicit artificial dissipation 
or of the upwind-type which are naturally dissipative. Both 
implicit and explicit temporal discretizations are discussed for 
the time-integration of the governing fluid flow equations. De- 
tails on the spatial and temporal adaption procedures are also 



Unstructured surface grid for F/A IK fighter configuration 


Fig ! 

given. The selected results that are presented demonstrate two- 
and three-dimensional applications for both steady and unsteady 
flows. Comparisons are also made with solutions obtained using 
a structured grid code and with experimental data to determine 
the accuracy of the unstructured grid methodology. 

Central - Difference Type Flow Solver 

The unsteady Fulcr equations in integral form are solved 
using a finite-volume algorithm that was developed for use on 
unstructured grids of triangles in 21) or ictrahedra in 3D. m 1-1 
The algorithm reduces conceptually to central differencing on a 
rectangular mesh and thus is referred to as a central difference 
type flow solver With this solver, artificial dissipation is added 
explicitly to prevent oscillations near shock waves and to damp 
high frequency uncoupled error modes. Specifically, an adap- 
tive blend of harmonic anti bihamionic operators is used, cor 
responding to second and fourth difference dissipation, respec- 
tively. The biharmonic operator provides a background dissipa- 
tion to damp high frequency errors and the harmonic operator 
prevents oscillations near shock waves. 

The Filler equations are integrated in time using a standard, 
explicit, four-stage, Rungc-Kutta time stepping scheme. In this 
scheme the convective operator is evaluated at each stage and, 
for computational efficiency, the dissipative operator is evalu- 
ated only at the first stage The scheme is second order accurate 
in time and includes the necessary terms to account for changes 
in cell volumes due to a moving or deforming mesh. Further- 
more, this explicit-scheme has a step size that is limited by the 
Couram-Friedncks-Ijcwy (CFL) condition corresponding to a 
CFL number of To accelerate convergence to steady-state, 
the CFF number may be increased by averaging implicitly the 
residual with values at neighboring grid points. These implicit 
equations are solved approximately using several Jacobi iteru 
lions. Convergence to steady state is further accelerated using 
enthalpy damping and local time stepping. The local time step- 
ping uses the maximum allowable step size at each grid point 
as determined by a local stability analysis. For unsteady appli- 
cations, however, a global time step is usually used because of 
the time -accuracy requirement. The maximum allowable global 


lime step may be increased to a value that is larger than that dic- 
tated by the CFL condition by using a time accurate version of 
the residual smoothing. Alternatively, temporal adaption may 
be used which involves a spatially varying lime step, as de- 
scribed in a subsequent section. 

Upwind-Type Flow Solver 

The unsteady Euler equations may be solved alternatively 
by using upwind differencing and either flux-vector or flux- 
difference splitting similar to upwind schemes developed for use 
on structured meshes. 10, ,5 ~ tg The present unstructured gTid al- 
gorithm is thus referred to as an upwind-type flow solver. The 
spatial discretization of this solver involves a so-called flux- 
split approach based on either the flux-vector splitting of van 
Ixrer 20 or the flux-difference splitting of Roe. 21 These flux-split 
discretizations account for the local wave- propagation charac- 
teristics of the flow and they capture shock waves sharply with 
at most one grid point within the shock structure. A further 
advantage is that these discretizations are naturally dissipative 
and consequently do not require additional artificial dissipation 
terms or the adjustment of free parameters to control the dissi- 
pation. 

The Filler equations are integrated in time using either an 
explicit Rungc-Kutta method (described in the previous section) 
or an implicit time-integration scheme involving a Gauss-Seidel 
relaxation procedure. )S The procedure is implemented by re- 
ordering the elements that make up the unstructured mesh from 
upstream to downstream. The solution is obtained by sweeping 
two times through the mesh as dictated by stability considera- 
tions. The first sweep is performed in the direction from up- 
stream to downstream and the second sweep is from downstream 
to upstream. For purely supersonic flows the second sweep is 
unnecessary. This relaxation scheme is unconditionally stable 
and thus allows the selection of the step size based on temporal 
accuracy of the problem being considered, rather than on the 
numerical stability of the algorithm. Consequently, very large 
time steps may be used for rapid convergence to steady state, 
and an appropriate step siz« may be selected for unsteady cases, 
independent of numerical stability issues. 
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Effects of adaptive mesh refinement on the total pressure loss contours for a 75° 
swept flat plate delta wing computed using the conical Euler equations at 
Af lVi = 1.4, o = 20°, and fl = 10°. 
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Spatial Adaption Procedure 

Spatial adaption is employed with the unstructured grid flow 
solvers to enrich the mesh locally in regions of high spatial flow 
gradients to resolve more accurately and efficiently the physics 
of the flow. 13 Equally attractive are coarsening techniques that 
remove elements from regions where relatively small changes 
in the flow variables occur. Both enrichment and coarsening 
procedures are currently being developed. 19 However, only the 
enrichment procedure is described herein. The enrichment pro- 
cedure uses an indicator to determine if an element in the mesh 
is to be refined or subdivided into smaller elements. Typically, 
the absolute change in density along an edge is used as an indi- 
cator for flows with shock waves and total pressure loss is used 
for flows with vortices. More recently a refinement indicator 
based on the material derivative of density 19 has been shown to 
be a superior indicator for unsteady flows. In general, the refine- 
ment indicator is compared with a preset tolerance to determine 
whether a given element should be refined. If the tolerance is 
exceeded, a new node is created at the midpoint of the edge 
and the element is divided. Each time the mesh is refined, an 
element may be divided in one of several different ways. The 
coordinates of the new node are determined by averaging the 
coordinates of the endpoints that make up the bisected edge. 
Special care must be taken, however, when an edge that is to 
be divided lies on a boundary of the grid, since the midpoint of 
the edge does not generally lie on the boundary. In this case, 
the location of the new node is determined generally by using 
a spline of the boundary coordinates. 

Temporal Adaption Procedure 

Temporal adaption is employed with the unstructured grid 
flow solvers, similar to spatial adaption, to resolve more accu- 
rately and efficiently the physics of the flow in time. 18 Tem- 
poral adaption can be thought of as time-accurate local time- 
stepping. Local time-stepping is typically used in a non-time- 
accurate manner to accelerate the convergence of the governing 
fluid flow equations to steady state. Since only steady-state is 
desired, it docs not matter that every point in the flow is at a dif- 
ferent time. This, of course, is not the case for a time-accurate 
problem, since each point in the flow for such a calculation must 
be on the same temporal level to maintain time-accuracy. The 
problem is that if all of the grid cells are marched at the same 
time step witli an explicit lime-marching scheme, the most re- 
strictive lime step must be used in order to maintain numerical 
stability. Temporal adaption is a method to march each cell at 
its own lime step, although ultimately the flow variables in all 
cells reach the same point in time. Temporal adaption can be 
viewed as similar to spatial adaption in that small lime steps 
should be taken only in localized areas governed by the flow 
physics and not in the entire flow field. Typically, small grid 
cells are integrated with small time steps and large grid cells 
arc integrated with large time steps. All of the cells reach the 
same lime level n + 1 to maintain time accuracy by using lo- 
cal time steps that are multiples of one another. The solution 
is integrated in a special sequence so that all values necessary 
for the calculations at an intermediate level are available at the 
proper times. For a particular cell to be integrated from time 
level n to time level u f 1, for example, the solution must also 
be known at its neighboring cells at time level n. If the value 
needed for the integration is unknown at a particular temporal 
mxle, it is determined from a linear interpolation between two 
known values. 


Deforming Mesh Algorithm 

For problems where the aircraft moves or deforms, the mesh 
must move so that it continuously conforms to the instantaneous 
shape or position of the vehicle. This is accomplished by 
using a spring network to model the original mesh such that 
each edge of the triangle or tetrahedron is represented by a 
spring. 12 The spring stiffness for a given edge is taken to be 
inversely proportional to the length of the edge. Grid points 
on the outer boundary of the mesh are held fixed and the 
instantaneous locations of the points on the inner boundary 
(aircraft) are given by the prescribed surface motion. At each 
time step, the static equilibrium equations in the x, y, and 
z directions, which result from a summation of forces, are 
solved iteratively at each interior node of the grid for the 
displacements. This is accomplished by using a predictor- 
corrector procedure, which first predicts the displacements of the 
nodes by extrapolation from grids at previous time levels and 
then corrects these displacements using several Jacobi iterations 
of the static equilibrium equations. The predictor-corrector 
procedure has been found to be more efficient than simply 
performing Jacobi iterations because far fewer iterations are 
required to achieve acceptable convergence. In practice it has 
been found that only one or two iterations are sufficient to 
accurately move the mesh. 

Results and Discussion 

Selected results from the unstructured-grid methods of Refs. 
10-19 are presented for two- and three-dimensional geometries 
for both steady and unsteady flows. Comparisons are made 
with solutions obtained using a structured grid code and with 
experimental data to determine the accuracy of the methodology. 

Two-Dimensional Euler Results 

To assess the two-dimensional central-difference-type Euler 
flow solver, calculations were performed for the NACA 0012 
airfoil. 11 These results were obtained using the unstructured grid 
shown in Fig, 4. The grid has 3300 nodes, 6466 triangles, and 
extends 20 chordlengths from the airfoil with a circular outer 
boundary. Also there are 1 10 points that lie on the airfoil sur- 
face. Generalized aerodynamic forces for the NACA 0012 air- 
foil oscillating in either plunge or piich-about-the quarter-chord 



Fig. 4 Partial view of unstructured grid of triangles 
about the NACA 0012 airfoil. 
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o Euler * harmonic 
4 r □ CFL3D - harmonic 



(a) due to plunge. (a) convergence histories. 
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(b) due to pitch. 

Fig. 5 Comparisons of generalized aerodynamic forces 
computed using CFL3D and the unstructured- grid 
central -difference -type Euler flow solver for the 
NACA 0012 airfoil at A/,*. = 0.8 and o 0 = 0°. 
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Fig. 6 Comparisons of generalized displacements 

computed using CFL3D and the unstructured-grid 
central difference-type Euler flow solver for the 
NACA 0012 airfoil at M^. = 0.8 and « p = 0°. 
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(b) steady pressure distribution. 

Fig. 7 Comparisons of steady-stale results for the 

NACA 0012 airfoil at Moo = 0.8 and a 0 = 1.25° 
computed using the upwind-type Euler flow 
solver with flux-vector splitting. 

are shown in Fig. 5. For this case the freestream Mach number 
was A/qc =0.8 and the angle of attack was a 0 = 0°. Com- 
parisons are given among results obtained using the so-called 
pulse transfer-function analysis, the harmonic analysis, and a 
harmonic analysis performed using a structured grid Navier- 
Stokcs code (CFL3D) run in an Euler mode. The three sets of 
results agree well with one another, for both plunge and pitch 
motions, thus verifying the accuracy of the unstructured grid 
method. 

Aeroelastic results were also obtained for the NACA 0012 
airfoil with two degrees-of- freedom (pitch and plunge) at Moo- 
0.8 and n Q = 0 0 . 11 Comparisons of second mode generalized 
displacements (q 2 ) are shown in Fig. 6 for three values of 
nondimensionai dynamic pressure ( Q ) that bracket the flutter 
point. The generalized displacements agree well with the struc- 
tured grid (CFL3D) solution which verifies the unstructured grid 
methodology for aeroelastic analysis. The flutter speed for this 
case, determined by interpolation of the dominant damping of 
these generalized displacements, also agrees to within 2% of the 
CFL3D value. 
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To test the more -recently developed upwind type Euler flow 
solver, steady How results were obtained for the NACA 0012 
airfoil at M ^ - 0.8 and u„ = 1.25°, using both implicit and ex 
plicil time marching. 1 ** The explicit rime marching results were 
obtained using a CFL number of 2.5 and the implicit time- 
marching results were obtained using a CFL number of 100,000. 
A comparison of the convergence histories is shown in Fig. 7(a) 
and the resulting steady pressure distribution is shown in Fig. 
7(b). The “error” in the solution was taken to be the L,»-norm 
of the density residual. As shown, the explicit solution is very 
slow to converge whereas the implicit solution is converged 
to four orders of magnitude in only approximately 500 steps. 
Also, the pressure distributions indicate that there is only one 
grid point within the shock structure, on cither the upper or 
lower surface of the airfoil, due to the sharp shock capturing 
ability of the flux splitting. Converged steady solutions arc thus 
obtained with the implicit algorithm with an order of magni- 
tude less CPU time than the explicit algorithm, and the shock 
waves are more sharply captured with the flux-split spatial dis- 
cretization than the central-difference* type discretization. These 
improvements in accuracy and efficiency arc also realized for 
unsteady applications. 

Conical Fuler Results 

Calculations were performed using the conical Euler version 
of the central-difference-type flow solver to investigate unsteady 
vortex-dominated flows about highly-swept delta wings. 14 This 
code includes the additional analysis of the free-to-roll case by 
the inclusion of the rigid-body equation of motion for simultane- 
ous time integration with the governing flow equations. Results 
were obtained for a 75° swept delta wing at a freestream Mach 
number of 1.2 and an angle of attack of 30°. A partial view of 
the grid which was used is shown in Fig. 8. The grid, which 
has a total of 4226 nodes and 8299 elements, indicates that the 
wing has thickness and sharp leading edges. 

Figure 9 shows the free-to-roll response of the wing which 
was initiated by using an initial angular velocity. In this cal- 
culation, for simplicity, the mesh was moved to conform to the 
instantaneous position of the wing by rotating rigidly according 
to the wing roll angle, rather than by using the deforming mesh 
algorithm. The results indicate that initially the oscillatory re- 
sponse diverges for small values of roll angle. As the angle 



Fig. 8 Partial view of unstructured grid about a 75° 
swept delta wing. 


increases to around 35°, the rate of divergence decreases due 
to stabilizing aerodynamics, and finally, the response reaches a 
maximum amplitude of motion at <f> = 38° corresponding to a 
limit cycle. These results are similar in_nature to those obtained 
by Arena and Nelson 22 in a low-speed experimental investiga- 
tion of wing rock. The wing-rock rime history from Ref. 22, 
shown in Fig. 10, was obtained for an 80° swept delta wing at 
30° angle of attack. Although the case considered in Fig. 9 is 
different from that of Ref. 22 (the data from Ref. 22 are for 
low speed flows whereas the conical Euler code is limited to 
supersonic freestream applications), the similarity between the 
two sets of results in Figs. 9 and 10 is noteworthy and gives 
credibility to the present calculations. 

Three-Dimensional Euler Results 

Unsteady flow results were obtained for a supersonic fighter 
configuration that was oscillating in a complete-vehicle bending 
mode to demonstrate a three-dimensional application of the 
central-difference-typc Euler solver. 12 The results were obtained 
using a grid which has 13,832 nodes and 70,125 tetrahedra. 
The surface triangularion of the aircraft is shown in Fig. 1 1(a) 
and the bending mode shape (exaggerated by a factor of five) 
is shown in Fig. 11(b). Instantaneous pressure distributions 
on the surface of the vehicle at the maximum (bend-up) and 
minimum (bend-down) amplitudes of oscillation are shown in 
Fig. 12. For this case the freestream Mach number was 
2.0, the reduced frequency (based on wing tip semi-chord) 
was 0.1, and two angles of attack of 0 and 12 degrees were 
considered. The results of Fig. 12 show the effects of angle 
of attack on unsteady pressures, and clearly demonstrate that 
the unstructured grid methodology can treat complex aircraft 
configurations undergoing structural deformation. 



Nondlmenslonal Time 


Fig. 9 Free-to-roll rime history for a 75° swept delta 
wing at ;l/ x = 1.2 and a - 30°. 



Fig. 10 Wing-rock time history for an 80° swept delta 
wing at 30° angle of attack (Ref. 22, reprinted 
with permission from Professor Robert C. Nelson, 
Notre Dame University). 
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To test the more recently-developed upwind-type Euler flow 
solver, calculations were performed for the ONERA-M6 wing. 16 
The M6 wing has a leading edge sweep angle of 30°, an 
aspect ratio of 3.8, and a taper ratio of 0.562. The airfoil 
section of the wing is the ONFRA “D” airfoil which is a 10% 
maximum thick ness -to chord ratio conventional section. The 



(a) original surface grid. 

Fig. 1 1 Surface grid for the 
configuration. 

• a = 0 ° 



results were obtained using a grid which has 42,410 nodes 
and 231,507 tetrahedra. Results were obtained for the M6 
wing at a freestream Mach number of 0.84 and 3.06° angle 
of attack. These conditions were chosen for comparison with 
the experimental pressure data of Ref. 23 as shown in Fig. 13. 
Hie results indicate that there is a weak supersonic-to-supersonic 



(b) assumed bending mode. 


Langley supersonic fighter 

• a = 12 ° 




Fig. 12 Effects of angle of attack on the instantaneous pressure coefficient contours on the 
I .angle y lighter configuration at the maximum (bend-up) and minimum 
(bend-down) amplitudes of deformation computed using the 
central-difference- type Filler flow solver at M tXj - 2.0 and k = 0.1. 
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shock wave in the inboard region, forward toward the leading- 
edge. The primary, supersonic- to subsonic shock which occurs 
in the midchord region coalesces with the first shock in the 
outboard direction toward the wing tip. Near the tip, the two 
shocks merge to form a single, strong, supersonic-to-subsonic 
shock wave. The Euler results are in fairly good agreement 
with the experimental pressure data, especially in predicting the 
strength and location of the shock waves, which tends to verify 
the upwind-type algorithm. The shocks are sharply captured 
with only one grid point within the shock structure, due to the 
flux splitting. 

Spatial Adaption Results 

To demonstrate the spatial adaption procedures, results are 
obtained for the NACA 0012 airfoil pitching harmonically about 
the quarter chord. iy The freestream Mach number was 0.755 and 
the mean angle of attack was 0.016°. The pitching amplitude 
was 2.51° and the reduced frequency (based on semi-chord) 
was 0.0814. Figure 14 shows the instantaneous adapted me s lies 
and Fig. 15 shows the corresponding instantaneous density 
contour lines (Ap = 0.02). The instantaneous meshes and 
density contour lines during the third cycle of motion were 
plotted at eight points in time. In each plot, the instantaneous 
pitch angle o(r) and the instantaneous angular position kr in 
the cycle are noted. The instantaneous meshes (Fig. 14) clearly 
indicate the enrichment in regions near the shock waves and near 
the stagnation points. They also show coarsened regions where 
previously enriched regions have relatively small flow gradients. 
The density contours during the cycle (Fig. 15) demonstrate 
the ability of die spatial adaption procedures to produce sharp 
transient shock waves. 
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Fig. 13 Comparisons of steady pressure distributions for 
the ONERA M6 wing computed using the 
upwind-type Euler flow solver with flux vector 
splitting at A/, x = 0.84 and n a = 3.06°. 


Temporal Adaption Results 

To demonstrate the efficiency of temporal adaption over 
global time-stepping for unsteady problems, results were ob- 
tained for the same NACA 0012 pitching airfoil case of the 
previous section. 18 Figure j 6 shows calculated results obtained 
using temporal adaption and global time -stepping as well as 
comparisons with the experimental pressure data of Ref, 24. In 
each pressure plot the instantaneous pitch angle o(r) and the 
angular position in the cycle kr are noted. During the first pan 
of the cycle there is a shock wave on the upper surface of the 
airfoil, and the flow over the lower surface is predominately 
subcritical. During the latter part of the cycle the flow about the 
upper surface is subcritical, and a shock forms along the lower 
surface. The pressure distributions indicate that the shock posi- 
tion oscillates over approximately 25% of the chord along each 
surface, and in general, that the two sets of calculated results 
compare well with each other This good agreement verifies the 
time -accuracy of the solution computed using temporal adap- 
tion, which was obtained at one -fourth of the CPU time that the 
global time-stepping solution required. Both sets of calculated 
results also agree well with the experimental data. 


Concluding Remarks 

The current status of unstructured-grid methods develop- 
ment in the Unsteady Aerodynamics Branch at NASA Lang- 
ley Research Center was described. These methods are being 
developed for unsteady aerodynamic and aeroelastic analyses. 
The paper highlighted the flow solvers that have been devel- 
op for the solution of the unsteady Euler equations and gave 
selected results which demonstrated various features of the ca- 
pability. The results demonstrated two- and three-dimensional 
applications for both steady and unsteady flows. Comparisons 
of two dimensional steady and unsteady results were made with 
solutions obtained using a structured grid code and with exper- 
imental data to determine the accuracy of the two dimensional 
flow solvers. Comparisons of three-dimensional steady results 
were also made with experimental data to determine the accu- 
racy of the three-dimensional flow solver. These comparisons 
showed good agreement which thus verifies the accuracy of the 
unstructured grid methods. 
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Fig. 16 Comparison of instantaneous pressure distributions for the NACA 0012 airfoil 
pitching harmonically at A/ x , = 0.755, = 0.016°, tt| = 2.51°, and k = 0.0814 

computed using the upwind type Euler flow solver with flux-vector splitting. 
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